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Higgs mediated flavor violating electromagnetic interactions, induced at the one-loop level by a 
nondiagonal Hfifj vertex, with fi and fj charged leptons or quarks, are studied within the context 
of a completely general effective Yukawa sector that comprises SUl{2) x (7i'(l)-invariant operators 
of up to dimension-six. Exact formulae for the one-loop 7/i/j and yyfifj couplings are presented 
and their related processes used to study the phenomena of Higgs mediated lepton flavor violation. 
The experimental limit on the /i — > 67 decay is used to derive a bound on the branching ratio of the 
/i 677 transition, which is 6 orders of magnitude stronger than the current experimental limit. 
Previous results on the t —> jj."/ and r fi-y^ decays are reproduced. The possibility of detecting 
signals of lepton flavor violation at 77 colliders is explored through the 77 — » klj reaction, putting 
special emphasis on the final state. Using the bound imposed on the Hrfi vertex by the current 
experimental data on the muon anomalous magnetic moment, it is found that about half a hundred 
events may be produced in the International Linear Collider. 
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I. INTRODUCTION 



The International Linear Collider (ILC) [1,] is an ambitious program of electron-positron collisions at the TeV scale, 
which will provide a clean environment to make studies beyond the capabilities of the Large Hadron Collider(LHC). In 
addition to complementing and adding precision to the LHC discoveries, the ILC will provide important information 
to our knowledge of the standard model (SM) and promises access to new physics that could eventually show up at 
the TeV scale. In particular, the operation of this collider in the 77 mode offers a unique opportunity to explore new 
physics effects through production mechanisms that are not accessible in leptonic or hadronic machines. This collider 
will reveal crucial information on those production mechanisms that are naturally suppressed in electron-positron 
collisions , as the involved cross sections can be significantly larger than the corresponding e+e~ ones. This is the case 
of processes that lead to flavor violation, which are quite suppressed within the standard model. In the quark sector 
of the SM, flavor changing neutral current transitions are very suppressed [3, whereas the analogous transitions 
in the lepton sector are absent at any order of perturbation theory. Although absent in the SM, the phenomena of 
lepton flavor violation (LFV) can arise in many of its well-motivated extensions. One new and interesting ingredient 
of most models beyond the SM is the presence of more complicated Yukawa sectors, which naturally favor nondiagonal 
interactions mediated by the physical Higgs bosons of the theory. Thus, extended Yukawa sectors are a good place to 
search for LFV transitions. In particular, one interesting production mechanism would be the one leading to the LFV 
T/ii final state. However, if this transition is mediated by a Higgs boson, one expects an uninteresting cross section 
for the e+e" Tfj, production mechanism at the tree level and beyond. The reason is that these types of processes 
that involve couplings among Higgs bosons and leptons are naturally scaled by the masses of the leptons. In this 
way, the cross section for the e"'"e~ Tfj, process is necessarily small, as it is proportional to the product of the three 
leptonic masses. Although first generated at the one-loop level, the 77 —^ mechanism offers a more interesting 
scenario. One advantage consists in the fact that the cross section depends only of the muon and tau masses, which 
are substantially larger than the electron mass. In addition, as already mentioned, in many cases the 77 cross sections 
can be much more larger than those associated with e~^e~ collisions. 

In this paper we are interested in investigating Higgs mediated fiavor violating electromagnetic transitions through 
the one-loop induced jfifj and "tjfifj couplings and their related processes, namely, the fi — > fjj and fi — > fj"f"f 
decays, as well as the 77 — *■ fifj reaction in the context of 77 colliders. Although we will present exact expressions 
for these processes, valid for quarks of leptons, we will focus on their implications on LFV. Although, as already 
mentioned, these types of processes are absent at any order of perturbation theory within the SM, they can be 
induced in many of its well-motivated extensions. In particular, we will exploit the potential of the 77 reaction 
at 77 colliders to search signals of LFV. We will focus on extended Yukawa sectors that are always present within the 
SM with additional S'[/L(2)-Higgs multiplets or in larger gauge groups. Some processes naturally associated to fiavor 
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violation could be significantly impacted by Yukawa sectors associated with multi-Higgs models, as it is expected 
that more complicated Higgs sectors tend to favor this class of new physics effects. We will assume that the jfifj and 
llfifj couplings are generated by a virtual scalar field with a mass of the order of the Fermi scale v « 246 GeV through 
the flavor violating H fifj vertex. However, instead of tackling the problem in a specific model, we will adopt a model 
independent approach by using the effective Lagrangian technique 4] , which is an appropriate scheme to study those 
processes that are suppressed or forbid in the SM. As it has been shown in Refs. [5, 6], it is not necessary to introduce 
new degrees of freedom in order to generate flavor violation at the level of classical action; the introduction of operators 
of dimension higher than four will be enough. We will see below that an effective Yukawa sector that incorporates 
S'[/i(2) X J7y (l)-invariants of up to dimension-six is enough to reproduce, in a model-independent manner, the main 
features that are common to extended Yukawa sectors, such as the presence of both flavor violation and CP violation. 
Although theories beyond the SM require more complicated Higgs sectors that include new physical scalars, we stress 
that our approach for studying flavor violating electromagnetic couplings induced by a relatively light scalar particle is 
sufficiently general to incorporate the most relevant aspects of extended theories, as in most cases, it is always possible 
to identify in an appropriate limit a SM-like Higgs boson whose couplings to pairs of W and Z bosons coincide with 
those given in the minimal SM. Besides its model independence, our framework has the advantage that it involves 
an equal or even smaller number of unknown parameters than those usually appearing in specific extended Yukawa 
sectors. As already mentioned, one important goal of this work is to calculate the amplitudes for the "ffifj and 'yjfifj 
interactions, which are generated at the one-loop level by the flavor violating H fifj vertex. We will present exact 
formulae, which means that the masses of all particles will be conserved, no approximations will be made. These 
general expressions can be used in searching signals of flavor violation both in the quark sector and in the Icpton 
sector. In a previous paper by some of us [6[, the implications of Higgs mediated LFV on the tau decays r — > /Z7 
and r — s- were investigated within the context of this model-independent description of extended Higgs sectors. 
In the present paper, we will concentrate in investigating the potential of extended Yukawa sectors to induce lepton 
flavor violation via the electromagnetic fi ej and /x —^ 677 transitions or through the 77 Tfj, reaction, which is 
the most promising one to detect possible signals of LFV at the ILC. 

The paper has been organized as follows. In Sec. |TT1 an effective Lagrangian for the Yukawa sector that induces 
the Hfifj vertex is presented. Section [Hi] is devoted to calculating the exact amplitudes for the 'ffifj and 'yjfifj 
couplings. In Sec. IIVI the branching ratios for the li — > Ijj and h —^ Ij^j decays are studied. In particular, the 
experimental limit on the fi ej decays will be used to impose a bound on the fi — > 677 transition. Section |V] is 
devoted to calculating and discussing the cross section for the 77 Tji reaction in the context of the ILC. Finally, in 
Sec. IVII the conclusions are presented. 



II. EFFECTIVE LAGRANGIAN DESCRIPTION OF HIGGS MEDIATED FLAVOR VIOLATION 

In the SM the Yukawa sector is both flavor- and CP-conserving, but these effects can be generated at the tree 
level if new scalar fields are introduced. One alternative, which does not contemplate the explicit introduction of 
new degrees of freedom, consists in incorporating into the classical action the virtual effects of the heavy degrees by 
introducing SUl{2) x [/y (l)-invariant operators of dimension higher than four 5, 6]. Indeed, it is only necessary to 
extend the Yukawa sector with dimension-six operators to induce the most general coupling of the Higgs boson to 
quarks and leptons. A Yukawa sector with these features has the following structure [5|, l6|]: 

rf// = -F/,(i.$?,)-^($tcf>)(Z,ci>z^) + H.c. 

-r,^(Q.$d,) - ^(ci>tcO)(Q,<i>d^.) + H.c. (1) 

-Y^iQ.^Uj) - ^(ci>t$)(Q,$ii^.) + H.c, 

where Ky, Li, Qi, $, li, di, and Ui stand for the usual components of the Yukawa matrix, the left-handed lepton 
doublet, the left-handed quark doublet, the Higgs doublet, the right-handed charged lepton singlet, and the right- 
handed quark singlets of down and up type, respectively. The aij numbers are the components of a 3 x 3 general 
matrix, which parametrize the details of the underlying physics, whereas A is the typical scale of these new physics 
effects. 

After spontaneous symmetry breaking, this extended Yukawa sector can be diagonalized as usual via the unitary 
matrices V^''''" and V^''''", which relate gauge states to mass eigenstates. In the unitary gauge, the diagonalized 
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Lagrangian can be written as follows: 

(l + I^h) (eMiE + DMdD + UMuU) 

h(i + -^h(3 + -^h\\{eQ}PrE + DVL'^PrD + U^VPrU + H.c.) , (2) 
V Amw \ 2m\\r JJ\ J 

where the (a = l,d,u) are the diagonal mass matrix, whereas E — {e,fl,f), D = (d, s, 5), and U = {u,c,t) are 
vectors in the flavor space. In addition, fl'^ are matrices defined in the flavor space through the relation 

f^'^ = ^(|)Vi:aX'. (3) 

To generate Higgs-mediated flavor violation at the level of classical action, it is assumed that neither Y''''^'^ nor a''*^'" 
are diagonalized by the V£ ^ rotation matrices, which should only diagonalize the sum y'^'^."_|_Q,^'''", As a consequence, 
mass and interactions terms would not be simultaneously diagonalized as occurs in the dimension-four theory. In 
addition, if fi"^ ^ 17°, the Higgs boson couples to fermions through both scalar and pseudoscalar components, which 
in turn could lead to CP violation in some processes. As a consequence, the flavor violating coupling H fifj, where / 
stands for a charged lepton or quark, has the most general renormalizable structure of scalar and pseudoscalar type 
given by 

- iVij = -i{uJijPji + uj*jPl), (4) 

where uiij = ^ij + ^ij^ Pr — ^^2^ ; ^^^d P^ — ^-3^. This is the vertex that will be used in the next section to 
calculate the amplitudes for the ^fifj and 'yyfifj couplings, which are induced at the one-loop level. 

To close this section, let us to emphasize that the above effective Lagrangian describes the most general coupling 
of renormalizable type of a scalar field to pairs of fermions, which reproduces the main features of most of extended 
Yukawa sectors, as the most general version of the two-Higgs doublet model (THDM-III) [7] and multi-Higgs models 
that comprise additional multiplets of SUl{2) x C/f(1) or scalar representations of larger gauge groups. Our approach 
also covers more exotic formulations of flavor violation, as the so-called familon models or theories that involve 
an Abelian flavor symmetry ,9] . In this way, our results will be applicable to a wide variety of models that predict 
scalar-mediated flavor violation. 



III. THE ONE-LOOP AMPLITUDES FOR THE 77/«/j AND 7/,/^ COUPLINGS 

In this section, we present the amplitudes associated with the on-shell 77/i/j and jfifj couplings, which can be 
used to predict the cross sections or decay widths of their related processes. No approximations will be made and 
exact formulae will be presented. In particular, the amplitude for the 77/i/j coupling will be given in terms of 
scalar products among the four-vectors involved. From these expressions, the conversion to the kinematical variables 
associated with the scattering 77 fifj process (Mandelstam's variables) or with the three-body fi fjj-f decays 
(phase space variables) can be easily performed. 



A. The yyfifj coupling 

The contribution of the flavor violating H fifj vertex to the ^^fifj coupling is given through the box and reducible 
diagrams shown in Figs. [1] and [21 each leading separately to a finite and gauge-invariant result. It should be noticed 
that the contribution of the H fifj vertex occurs at the one-loop level in graphs of Fig. [Tl whereas it contributes at 
the tree-level in the diagrams of Fig. [51 This last contribution, which is determined by the one-loop i?77 vertex, can 
eventually have a determinant role in the cross section due to an resonant effect of the Higgs boson. In the case of 
those processes involving light fermions, such as the LFV 77 Tfj, reaction, the contributions arising from diagrams 
of Fig. [1] are marginal due to the presence (at first order in the fir^ parameter) of the SM coupling Htt, which is 
proportional to rriT/v. However, the contributions given by these diagrams can be competitive with that induced 
by the diagrams of Fig. [21 as this suppression effect is not present in flavor violating processes involving the quark 
top [lo| . We have adopted the convention of taking all momenta incoming: ki + k2 + Pi + Pj — 0, where fci 2 and pij 
stand for the four-vectors of photons and fermions, respectively. The rest of our notation is indicated in Figs. [T] and 

m 

The analytical structure of the amplitude is dictated by electromagnetic gauge invariance and Bose statistics. In 
addition, since the structure of the H fifj vertex is of renormalizable type, the amplitude must be free of ultraviolet 
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divergences. We first analyze the contribution given by the set of diagrams shown in Fig. [T] The corresponding 
amphtude can be written as foUows: 



■'^'Box+rod 



2 ^fJ-i^ 



^ A, 



n— 1 " 



6 



^ Ac 



(5) 



where the subindexes B, T, and 5 stand for box, triangle, and self-energy diagrams, respectively. The diverse 
magnitudes appearing in these expressions are given by 



-'t-i 

^3 



{ujkj Pr + ujI^ Pl) - i>j + rrik) Y + + ]/>i + nik) Y + ]/>i + mk){uJik Pr + uj*^ Pl), 
(iJkj Pr + ujIj Pl) - ^- + mk) 7^ {^ + h+ k + "^fc) l"" (1^ + A + "^fe)(w^fc Pr. + '^Ih Pl), 
(ujkj Pr + ujIj Pl) (/^ - ^- + nik) Y {^ + h+k+ ^k) (^^.fc Pr + Pl) (.h + A + '^^) 7^, 
(cjfcj Pr + w^j- Pl) - ^- + mk) 7^ + /^■2 + A + mfe) (t^ifc Pr + oj*^ Pl) {^2 + ^ + rui) 7^ 
7" i^i +k + m,) (LUkj Pr + Pl) {$ + ^1+1/), + mu) 7^ ^k){^^k Pr + <fe Pl), 

7^ (1^2 + A + mj) (tjfcj Pr + ljIj Pl) {1/^ + ^2 + i>i + mk) 7" {^ + i>i+ mk)iuJik Pr + uj*k Pl), 
7^ (-^- - /^i + mj) Y {i>^ + Pfl, + tJ^j- Pi) (1^ + A + mfe)(a;,fc + Pl), 

-i" {-i>^ -1^2+ TOj) 7^ (A + "^^) ^i?^ + <j Pl) + ^ + mk){uJ^k Pr + Pl), 
7^ {-^j mj) [ujkj Pr + i^lj PL){^-h- h + '^k) (^^-t + ^^fe ("^1 " + "^0 t'', 

7^^ - + mj) Pr + t^^^- Pl) {^-^2- ^ + mu) {uj^u Pr + u:*^ Pl) {-h - h + "^') 
Pr + t^;!^- Pl) (/^ - ^- + mk) {uJ^k Pr + uj*k Pl) {-^ + m,) -1" (/^i + + m,) 7^, 
Pi? + oj*kj Pl) - ^- + mk) {uj^k Pr + uj*^ Pl) {-^ + m,) 7^ {^2 + A + ^i) l"" ^ 



(6) 



and 



[(fc - ml\\{k + A:„ - ml\\(k-rV^f " - "4], 

r [(/c-pj)^-?'^fe][(fc + fcn+J5»)^-w2][(fc„+p,)2-mf][P_m|^], n = l,2 

\ [(A:„-2+P«)^-?'^?][(A: + /c„-2+ft)^-m2][(fc + p,)^-"i^][^^-"ili], n = 3,4 ' 

f [(Pj+A:„)2-m2][p2_TO2]p^p^)2_^2j[p_^2^]^ n = l,2 

< [(Pj+A:„_2)^-m2][(fc-pj-fc„_2)^-m2][(pj +fc„_2)^-mf][A:2-m|^], n = 3,4 . (7) 

[ [(A;-p,)2-m2][p2_^2][(;.^_^^p^)2_^2][fc2_^2^]^ ^ ^ 5^ g 

Gauge invariance means that the amplitude must satisfy the transversality conditions, 

fclM-A^rox+..ed(fcl-fc2)-0, (8) 
fc2.A^rox+rcd(^1.^2)=0, (9) 

whereas Bose statics requires a symmetric amplitude under the interchanges k\ k2 and /i <-!■ z/: 

MZ.+Uki,k2) = AlBl+.cd(^2, A:i). (10) 

To solve the above integrals, we have used the Passarino-Veltman Lorentz tensorial decomposition pH implemented 
in the FeynCalc computer program Although, in order to make less cumbersome the expressions, we present 

our results in terms of the form factors that define this Lorentz decomposition, we have verified that they are free of 
ultraviolet divergences by performing the reduction of these form factors to Bq, Cq, and Dq scalar functions. Although 
the box diagrams depend on Bq functions, the divergences cancel among themselves dealing to a finite amplitude. 
This is in contrast with the sets of triangle and self-energy graphs, which generate divergent pieces that disappear 
only after adding both sets of amplitudes. This cancellation of ultraviolet divergences is quite intricate due to the 
fact that the amplitudes depend on the Bq functions in a very complicated way. It is worth mentioning that a similar 
problem is present in the amplitude for the one-loop t egg decay of the quark top within the context of the SM. 
The amplitude that presents this problem is induced by the contribution of the pseudo-Goldstone boson associated 
with the W gauge boson, as it was studied in Ref. 0]. The authors of this reference verified numerically that their 
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results are free of divergences. In our case, we carried out an analytical analysis in order to convince ourselves that 
these divergences are absent indeed. To perform this analysis, we replace each of the Bq functions by its divergent 
part, which is common to all them, namely, we put Bq A — — 1/(Z? — 4) in our amplitude in order to express the 
apparently divergent part in a product of the way FR'^'^ A, with F a scalar function that depends in complicated way 
on kinematics variables and i?^"^ a nontrivial Lorentz tensor structure. After some nontrivial algebraic manipulations, 
we achieve the required result F = 0, which shows that in fact the amplitude is free of ultraviolet divergences. As 
to gauge invariance, it is achieved only after adding the amplitudes of all diagrams of Fig. [TJ On the other hand, 
the contribution from the set of diagrams shown in Fig. [2] is free of ultraviolet divergences and is gauge-invariant, as 
it is generated essentially by the SM one-loop Hjj vertex, which, as is well-known [l^l, is characterized by a finite 
and gauge-invariant amplitude. After these considerations, the complete amplitude for the 77/i/j coupling can be 
conveniently written in the following way: 

where TBox-i-rcd ^'^'^ ^red(H) represent the contributions coming from the sets of diagrams given in Figs. [T] and [H 
respectively. We have normalized the amplitudes to the W gauge boson mass mw- The rg^^^j.j,j tensor amplitude 
can be organized in terms of 11 gauge structures as follows: 

r^ox+red = E (Tt: Pr + Pl), (12) 

where gauge structures means that they satisfy the transversality conditions 

ki.Tt:, = 0, (13) 
k2uTt:^ = 0. (14) 

We have arranged our results in an order that reflects in a manifest way both gauge invariance and Bose symmetry. 
Accordingly, each of the T^^^^ quantities are made of products that contain a gauge structure and a scalar form factor, 
which respects both symmetries as a whole. These tensor amplitudes are given by: 



'^-^ ki-k2 

_ (Pj fci ■ fca - fc^ fci ■ pj){p^ fci ■ fcz - k'( k2 ■ Pj) 

" ^^'^ {mw ki ■ k2f 

^ (7^ kr-k2-UiUMk^-k2-k\ k2-Pj) fci ■ fc2 - k- ^2){P^ k,-k2- k^ fci • p,) 

mw (fci • fc2)2 + mw (ki ■ k2f 

rrt^'^ -IT? ;/ , 17 APj ki ■ k2 - kl^ ki ■ pj){p'^ ki ■ k2 - k'( k2 ■ Pj) 

(fci • k2y 



Dr. 7? ^L,R 



7'^ (pJ fci • fc2 - fcr k2 -Pj) ^ h 7" (p; fci • fc2 - fc^' fci • V, ) 

—2 — z — z 1" P'^-^L.R —2 — z — r 



m^w ■ 


k2 




7^ ^2 Y kx-k2- 


- h h 1" 


fcj 


mw fci 


■k2 




hrr ki-k2' 


- h 7^ h 


fcr 



mL fci • fc2 



_^ 1-^2 1 n-1 • "-2 - ^l ^2 7 '-2 , J. Yhrf^i->^2-h^ii''K 

'^.R • m,.r i-1 . /c-o ' mw fci • fc2 

rr,f,u r? n r / a-i • n,2 - ^1 ,■ ^2-1 , /^2 7^ 7^ ^1 ' ^2 " /^2 7^ /^l ^2 

^8lr -^i^L.R ; ; <- , 

* mw fci • fc2 mw fci • fc2 

_ h hriPjki-k2~k^ki-p,) ^2^i7^(p,"fci-fc2-fcrfc2-Pj-) 

m3,fcl-fc2 m3,fcl.fc2 

_ ^I7^^2(pjfcl-fc2-fcrfc2-p,) |^27'^^l(P,^fcl-fc2-fc2^fcl-p,-) 

Jio..«-^i7.,« m3,fci.fc2 m3^fci.fc2 

^llL.R -^19l.r —2 ^ ^20l.r —2 ■ 
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A,(*:a),/l„(ti 




FIG. 1: Box and reducible diagrams contributing to the 77/i/j coupling. 

(i-.i ■ — 




-4„(fa) 



A^(h) 

A^ih) 




FIG. 2: Contribution of the standard model loop-induced ^77 vertex to the 'yyfifj coupling. 

The fact that each T^^ ^ term respects the Bose symmetry means that, for example, i^ig^ ^ transforms into ^20^ l 
under the interchanges fci <-> k2 and ^ v, but i^i^ remains invariant under the same interchanges. On the other 
hand, the rfgd(H) amplitude is given by 



cd(H) 



2 fci • /c2 - mfj + i WLH Th 



(16) 



The form factors i^„^ ^ are presented in the Appendix. 



B. The 7/,/j coupling 

The contribution of the Hfifj vertex to the ^fifj coupling is given by the diagrams shown in Fig. [31 The amplitude is 
entirely governed by electromagnetic gauge invariance. In principle, the amplitude could depend on Lorentz structures 
of the following type: monopole (7^), anapole (g^7^ — ^9/^): magnetic dipole (cr^^g'^), and electric dipole {'^b(^^ii,q^)- 
However, since this is a flavor violating transition, there is no contribution to the monopole. Also, the anapolar 
structure is absent in this case because the three particles are on-shell. So, the amplitude must be made of two 
pieces, one CP-conserving associated with a magnetic dipole of transition and the other being CP-violating, which 
is associated with an electric dipole of transition. The latter amplitude only can arise if an imaginary part of the 
parameters is assumed. Since the anomalous Hfifj vertex has a renormalizable structure, these amplitudes must be 
free of ultraviolet divergences by renormalization theory. As in the previous case, we will take all momenta incoming, 
ie., Pi + Pj + q = 0, with q the photon four-vector. Then, the vertex associated with the jfifj coupling is given by 



r,. 




(17) 



7 






FIG. 3; Diagrams contributing to the 7/i/j vertex. 



where 



(18) 
(19) 
(20) 



A3 = [k' ~mjj][ik + p,f -ml]. 



(21) 
(22) 
(23) 



Once the integrals are solved, one can see that the dependence on 7^ cancels exactly after using the Gordon identity, 
which allows us to eliminate the dependence of the amplitude on the four- vector p^^. The result can be written in a 
symmetric form as follows: 



where 



la 



j\(^Xij Xj ^ Xk) — {^i^ki -^j^kj^ {"^i T -^j T k) : 

A{xi,Xj,Xk) = {xiSki + Xj6kj)F{x.i,Xj,Xk), 



(24) 



(25) 
(26) 



F{xi , Xj , X}- ) — 



- I - Xk{xi ^ Xj + Xk)2m%Co + 



xf - xj 



(So(l)-Bo(2) 



1-4 

2 2 . 

j J \ .7 



^(i3o(3)-So(2)) -^(i?o(3)-i3o(l)) 

X j \ / Xq^ \ / 



(27) 



2(.Xj + Xj) 

' '^^(Bo(2)-i3o(3) 



1 I -^i 



2 2 1 



1 - Xk{xi - Xj + Xk)2m'jfCo + 

^(i3o(l)-i?o(3) 



2(1 - 2x1) - 2xk{xi - Xj) - XiXj 



2 2 



(so(l)-Bo(2) 



(28) 



In the above expressions, fi^ = Re{flij), ilfj — Im{ilij), and svK = sin^vi^. In addition, the dimension- 
less variables Xa = ma/mn were introduced. Also, Co — Co{'mf,m'j,0,m1,m'jj,m1), Bo{3) — Bo{0,m'jj,m'l), 
Bo{l) ~ Bo{nif jTJi'jjjnil), and Bq{2) — Bo{m'j,m'jj,m1) are Passarino-Veltman scalar functions. 
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IV. THE k Ijl AND h -> Ij-n DECAYS 



In this section, we analyze the branching ratios for the li Ij"/ and li —^ Ij"/"/ decays and confront our prediction 
with current experimental limits. We will concentrate in the n ^ ej and /i — s- 677 transitions, as those involving the r 
lepton were already studied in Ref. @, in which the limits Br{T fij) < 7.5 x 10"^" and Br{T A'77) < 2.3 x 10~^^ 
were derived. 



A. Experimental constraints from k Ij'y and /i eee 

From Eqs. p4l [27l [28]) . it is straightforward to construct the branching ratio for the fi fjj decay. It can be 
written as: 

where Ti is the total decay width of the fermion fi. This is an exact result, which can be used for studying Higgs 
mediated flavor violation in both the quark and the lepton sectors, but in this paper we will focus on LFV transitions. 
In particular, it is interesting to use the current experimental limits on the LFV li — *■ Ij^ transitions to constrain the 
riij parameters. We will focus on the LFV transitions of the muon. It is important to comment at this point that the 
complex component of fi^e also contributes to the electric dipole moment (de) of the electron. It can be shown that 
the current limit on induces a very strong bound on |14| , so from now on this contribution will be neglected. 
The limit on the branching ratio of the /i — > 67 decay reported by the Particle Data Group [ISj is 



Br^^pifi ^ 67) < 1.2 X 10-". (30) 

Then, working out at first order in the fl^^ = fJ^e parameter, one obtains the following constraint 

n^e<7x IQ-^. (31) 

This bound was obtained assuming a value for the Higgs mass of about the lower limit imposed by LEP on the Higgs 
mass 1(5;]. In the next subsection, we will use this constraint to impose a bound on the /i —^ 677 decay. 

It is interesting to compare this bound with the one that can be derived from the experimental limit on the branching 
ratio of the tree-body decay — > eee, which is strongly restricted by the experiment. In the context of our effective 
theory, this decay can occur through a virtual Higgs boson as it is shown in Fig. 2) The amplitude is governed by the 
LFV Hfie vertex and the SM one Hee, which is very suppressed. The corresponding branching ratio can be written 
as follows: 

i^K. - eee) = 7^) /(,,.), (32) 



where 



{x - yY + yz 





In the above expressions, y — rnj^/mj^ and z — r|^/m^, where F^ and Tjj are the muon and Higgs total decay 
widths, respectively. To bound the fJ^e parameter, we use the limit for Br{^ —^ eee) reported by the Particle Data 
Group [IBl, which is given by 

BvExpif^ eee) < 10'^^ (34) 

In Fig.[5]the behavior of lil^eP as a function the Higgs mass is displayed. It can be appreciated from this figure that 
\^fie\ < 0.3, at the best. This bound is almost 2 orders of magnitude less stringent than that derived above from 
the experimental limit on the — > e7 decay. The important difference between these bounds can be understood by 
noting that while the /i — *■ eee decay is governed by the LFV Hfie vertex and the SM one Hee, the ^ e^ transition 
is determined by the LFV Hjie vertex and the SM one if/^/i. Besides a phase space factor, the branching ratio for 
the /i — !■ eee decay is suppressed by an additional (me/m^)^ w 2.5 x 10~^ factor with respect to the one associated 
to the two-body ^ ^ ej decay. 
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B. The li Ij'yy decay 



The kinematical of the U —> ljj"f decay is characterized by three dimensionless variables x, y, and z (x + y + z ~ 2), 
which are related to the involved four-vector scalar products as follows: 



ki 


■Pi = 


k2 


■Pt = 


Pi 


Pj = 


ki 


k2 = 


ki 


Pj = 


k2 


Pj = 



2 ' 

2 



y m 

2~' 
z mf 
2~' 

{x + y + £,-l)m1 
2 ' 
[l-j-y) mf 
2 

(1 - C-a;)TOf 



(35) 



where ^ = rnj/mf. 

The decay width for the three-body li — > lj"f"f transition is determined by the squared of the amplitude given in 
Eq. . The corresponding branching ratio is given by 

i+? y+ 

Br{h - /,77) = ^ ^ J dxjdy \M\\ (36) 

2VI y- 

where A4 = A^'^'^e* (fci, Ai)e* (^2, A2), with e*(fci,Ai) and e*(fc2,A2) the polarization vectors of the photons. The 
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FIG. 6: Behavior of Br{iJ, —> 677) as a function of the Higgs mass. 



integration limits are given by 



1 r 

2 

1 r 



(2 - x) + v/a;2 -4^ 



y- 



(2 - x) - -4^ 



(37) 
(38) 



The T /i77 decay has already been studied in Ref. [(|. Here, we are interested in obtaining a bound on the /i ejj 
transition using the limit on the f2^e parameter obtained above. The behavior of the corresponding branching ratio 
is shown in Fig. [6] as a function of the Higgs mass. It can be appreciated from this figure that this branching ratio 
ranges from approximately 10~^^ to 10~^^. Conservatively, we can assume the following bound: 



i?r(/i — > 677) < 10 



-16 



(39) 



which would be the highest branching ratio allowed by the lower bound imposed by LEP [16| on the Higgs mass. 
This means that the /i — > 677 transition may be undetectable if it is mediated by a Higgs boson. This bound for 
Br{fj, —^ 677) is almost 6 orders of magnitude more stringent than the experimental limit reported by the Particle 
Data Group [l^, which is 7.2 x 10"^^. In obtaining the above constraint, only the contribution of diagrams given 
in Fig. [2] were considered, as the contribution of the set of diagrams shown in Fig. [T] is insignificant in this process 
involving very light fermions. 



V. THE CROSS SECTION FOR THE 77 ^ PROCESS 



In this section, we analyze the possibility of detecting signals of LFV at the ILC through the 77 — » T/i reaction. 
The constraint on the Hfxe coupling derived in the previous section from the experimental limit on the fi ej decay 
is so restrictive that no signals can be detected via the 77 fxe process. On the other hand, to predict the cross 
section for the 77 — > r/i reaction, we will need to assume some value for the VIt^ parameter. As in the case of the 
muon decays studied above, we will assume that the imaginary part of O,-^ is very suppressed compared with its real 
part, so from now on flrfi stands for the real part of this parameter. There are various experimental limits [l5| which 
can be used to bound fir/^j but the best constraint arises from the anomalous magnetic moment of the muon, which 
leads to the following bound [T3 |: 

nl^ < 1.1 X 10-3, (40) 

for 115 GeV < niH < 200 GeV. 

The kinematic of the 77 — > r/i process is determined by Mandelstam's variables, which are given by 
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■k2 = 






k2 




ki 


Pr = 


k2 


■Pt = 



,2 

2 ■ 



2 

^2 



2 

The unpolarized cross section for this process can be written as 



(41) 



a(77 ry) - — / dt lAlp, (42) 

16 TT^ Jtmi„ 

where Al is the invariant amphtudc which is given by Eq. pip . In the above expression, the integration hmits are 
given by 

ima.(t™»„) = -2 I 1--^^^-^) (lT\/r^), (43) 

where 

^-( .^"^r^ )^ (44) 

On the other hand, the convoluted cross section, i.e., the cross section for the complete process e+e^ ^ 77 — > r/^, 
can be written as 

ct(s)=/ dz— !^ct(77-> r^), (45) 

(nV+»»^)/\/s 

where s = z'^s, with ^/s(\/l) the center of mass energy of the e"*'e~(77) collision, and — is the luminosity of the 
photons, defined as 

dz 



= 2z / -f,,e{y)fy/e{z'/y), (46) 



where the energy spectrum for the backscattered photon is given by [17| 

/7/e y) = T^TT 1 - y + -j n T + ~Tn \2 ' (47 

£»(x) L l-y x(l-y) X (l-y) -l 

with 

^(X) = (1 - - - 4) log(l + X) + ^ + - - ^2 - (48) 

\ X 2 X 2(1 + x)'^ 

In this expression, x — (4£'ot^o)/™e' where and i?o are the mass and energy of the electron, respectively; ujo is the 
laser-photon energy; and y represents the fraction of the energy of the incident electron carried by the backscattered 
photon. The optimum values for the Umax and x parameters are ymax ~ 0.83 and x = 2(1 + \/2)- 

A. Discussion 



We now turn to discuss our results within the experimental context of the ILC. Although this collider is intended to 
operate initially at a center of mass energy of 500 GeV with a luminosity of f-e+e- = 500 fb~^ it is contemplated 
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m„ [GeV] s"2 [GeV] 



FIG. 7: Contribution of diagrams of Fig. \T\to the cross section a{e^e jj ^ fir) as a function of the Higgs mass for 
v's = 500 GeV (left) and as a function of ^/s for niH ~ 115 GeV (right). 

to increase the energy up to about 1 TeV in a subsequent stage. In the foUowing, we wiU to carry out our numerical 
analysis in the interval 500 GeV< ^/s < 3000 GeV. We will discuss separately the contributions to this process 
coming from diagrams displayed in Fig. [1] and Fig. [2l To analyze numerically our results, we will work at first order 
in the anomalous parameter J^t/j- We will use the value determined by the upper limit imposed by the anomalous 
magnetic moment of the muon, namely, fi^^ < 1.1 x 10"'^ [l3 |. It turns out to be that those contributions coming from 
diagrams shown in Fig. [T] are proportional to a sum of two amplitudes, each characterized by the factors (mr/v)(flT-^) 
and {m^/v){^Tti)i which determine the intensity of the SM coupling Hlili and the LFV one Hlilj. Such amplitudes 
represent the two possible situations that can occur when the r or /i particles circulate in the loops. Both such 
possibilities are implicit in our general and exact calculation represented by the amplitude given by Eq. (|lip . but we 
will neglect the contribution associated with a virtual muon, as it is quite suppressed with respect to the corresponding 
tau contribution. 

In Fig. [7] the contribution from diagrams displayed in Fig. [T] is shown as a function of the Higgs mass and the 
center of mass energy of the collider, whereas in Fig. [5] the contribution arising from diagrams of Fig. [2] are shown. 
From these figures, it can be appreciated that the contribution from diagrams of Fig. [1] are suppressed by about 7 
orders of magnitude compared with the contribution given by diagrams of Fig. [21 For = 500 GeV, the respective 
contribution to the cross section can reach a value of 10~^ fh for mn ~ 140 GeV, which shows clearly the marginal 

role played by diagrams of Fig. [Tl However, this situation may be different in those processes involving heavier 

fermions, as the 77 tc reaction [lO|, with t and c the top and charm quarks, respectively. Finally, the number of 
events as a function of the Higgs mass and the center of mass energy of the collider is presented in Fig. O From this 
figure, it can be appreciated that up to 50 events can be produced with a luminosity of £e+e- = 500/6^^. 

VI. CONCLUSIONS 

In the minimal standard model the couplings of the Higgs boson to the remaining massive particles are thoroughly 
determined, which should be considered an outstanding feature of the model. However, there are well-motivated 
extensions of the model that require the presence of extended Higgs sectors. The extended theory is of course 
less predictive due to the proliferation of free parameters. The bonus is the appearance of interesting new physics 
effects such as lepton flavor violating transitions that are mediated by the physical scalars of the theory. In this 
paper, the potential of extended Yukawa sectors in generating flavor violation both in the lepton and the quark 
sectors was explored in a model independent manner using the effective Lagrangian approach. A Yukawa sector 
that includes SUlI'^) x C/y (l)-invariant operators of up to dimension-six, which reproduces the main features of 
extended Yukawa sectors, was considered. The one-loop contribution of a flavor violating Hfifj coupling, where 
fi and fj stand for charged leptons or quarks, to the on-shell 'yfifj and Jjfifj electromagnetic interactions was 
calculated. Exact formulae for the amplitudes associated with these couplings were presented and used to investigate 
electromagnetic lepton flavor violating transitions. The experimental limit on the branching ratio for the /i — > 67 
decay was used to impose the bound \ ft^e\ < 7 x 10~^ on the if/ie coupling. This bound allowed us to derive the limit 
Br{ii 677) < 10^^^, which is about 6 orders of magnitude more stringent than the corresponding experimental 
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FIG. 8: Contribution of diagrams of Fig. [2] to the cross section cj{e^e jj ^t) as a function of the Higgs mass for 
^/s = 500 GeV (left) and as a function of ^/s for iriH ~ 115 GeV (right). 
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FIG. 9: The number of events as a function of the Higgs mass for ^/s — 500 GeV (left) and as a function of ^/s for mn = 115 
GeV (right).. 



limit. The possibility of detecting signals of lepton flavor violation at the ILC through the reaction e+e^ 
was explored. It was found that with the projected luminosity of this machine up to half a hundred of fM^ events 
could be produced. Although in some processes this number of events could be insufficient to detect the signal, it 
is important to have in mind that the fi^ signal is a clean environment and is spectacular in the sense that the 
process is free of background, having thus the chance of being observed. 
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APPENDIX: FORM FACTORS OF THE 77/i/j COUPLING 



In this appendix, we present the expressions for the form factors -Fn^ . We present results only for the right- 
handed form factors -Fki,, as it can be passed to the left-handed ones via the interchanges UkiL^mk ^ ^ki'^mk ^^'^ 
^^lii^mk ^ '^kii-^mk- There are a total of 20 form factors of this class, but they are linked by the Bose symmetry, so 
one needs to list only 11 of them, the rest being obtained through the interchange ki ^ fc2. Also, the arguments of 
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TABLE L Arguments for the Passarino-Veltman scalar functions. The expression has been defined as: h{ki, k2) — 2 (fci • k2 
{ki + k2)-pj). 



the diverse Passarino-Veltman coefficients are listed in Tables HI HIl IIIII v IIVI These form factors are given by: 



F2„ 



= Aujikmw (mkUJkj (-Doo(l) + -Doo(2)) - (£'002(1) + £'002(2))) , 

^UHk^v^^ (m, iu,, (£22(1) + £22(2)) - m, u;*,^ (1^222(1) + £222(2))) + ^^^^ [^j {ml - ml) u;*,^ 

X {ki + k2) ■ Pj Boih) +4{ki-k2 + (ki + k2) ■ Pj) {ml + 2 (/ci • fcj + (fci + fcs) • pj)) ((mfe ivkj (Co(4) 

+ Ci(6) + C2(6)) + Lul^ (Co (4) + 2Ci(6) + Cii(4) + 2 (Ci2(4) + ^2(6)) + ^22(4)) ) fci • Pj 

+ {mkiukj (Co (3) + Ci(5) + C2(5)) + m, (Co (3) + 2Ci(5) + Cii(3) + 2 (Ci2(3) + C2(5)) + C22(3))) 

X k2 ■ Pj^ + £0(^4) (fci + fc2) ■ Pj (mj {{ml + ml- m^) + 2 mj m^ ujkj) + 2 {mj uj^^ + 2 Wfc Lo^j) 

X (fci • fc2 + (fci + fc2) - Pj)) |, 

F3« = -Auol^u^umlv (i?oo(l) + i?ooi(l) + ^001(2) +i?003(l)), 



where 



g(fci, fc2) = fci -pj k2 -pj (fci • fc2 + (fci + fc2) -Pj) {ml + 2 (fci • fc2 + (fci + fc2) ■ Pj)) , 
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TABLE 11: Arguments for the Passarino-Veltman tensorial coefficients with an index. 
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TABLE 111: Arguments for the Passarino-Veltman tensorial coefficients with two indices. 
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TABLE IV: Arguments for the Passarino-Veltman tensorial coefficients with three indices. 
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- w^^. (Co (2) + Ci(4) + C2(4)) - mk cukj (Co (2) + Ci(4) + C2(4)) + c^*^. (Ci(2) + Cii(2) + 2 Ci2(2) 
+ C2(2) + C22(2))) - mi LO*k (mk to^j (Co(2) + Ci(4) + C2(4)) + mj tOkj (Co(2) + Ci(2) + Ci(4) + Cii(2) 

+ 2C„(2) + a(2) + CM + C.(2))))] + ,,^.,^(,„;^,,,.,j(,„._,„;^,,^,,,) ["'. h 

X ((m^ - ™j - '^ifc - 2mimkUJ*k) + [mjjmiUJik - mi (mj+ml) a;*^ - 2m'jmkUJik) i^kj) -Bo(2) 
+ ruj [m^j - ml) {mj uj^j uJik - mi a;*^ (Jkj) -Bo(ii) - 2 (^2 m, a;*^ + (m| + 2 - 2 m|^) w^fe 



■ mj {mi LOik +2mk tOik) i^kj^ k2 ■ Pj Bo{2) 



ki- k2 + (fci + ^2) • Pj 
2g{ki,k2) 



ojlj UJik {ki + fe) • ft ( {mjj - ml) 



X Bo{ii) + Bo{u) (to| + ml-m%+2{ki-k2 + {ki + k2) ■ Pj 
2 ujlj ujtk 



k2 ■ Pj 



Cl2(4)+C2(6)+C22(4)+2 (Dl22(l)+Dl22(2)+D223(l)) fe'ft >, 
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19h 



"21 



SinikCOkj (£'2(l)+£>2(2)) 



4 (1)22(1) + 1)22(2))) 



mj - ml) 



X (All + k2) ■ Pj Bo{ii) + 4 (All • fc2 + {ki + ^2) • Pj) (m| + 2 (fci • A;2 + (fci + A;2) • pj)) ((m^ Wfe^ Co(4) 
+ mjijlj ((7o(4) + Ci(6) + C2(6)))fci -p^ + (mfeW^,- (Co(3) + Ci(5) + ^2(5)) + m,- w^^. (Co(3) + 2Ci(5) 
+ Cn(3) + 2 ((7i2(3) + C2(5)) + (^22(3))) ^2 • p,) + Bo{h) (fci + ^2) • Pj (mj ((m| + - m|,) ujI 



+ 2 rrife w^j) + 2 (m^ loIj + 2 to^ w^j) (fci • ^2 + (A:i + ^2) ■ Pj) 



Phr = 



i^kj ^ik IT^W 



2 (Boiti) - Bo{t2) - m\ Co(3) + ml Co(3) - m] (Co(3) + 2 (Ci(5) + ^2(6))) - 2C2(5) kx ■ k^) 



4 fci • Pj 

- 4 (Co(3) + Ci(5) + Ci(io) + C2(5) - 2£)oo(l) + 2i)i(l) /ci • ^2) fci • - 4^2(5) • Pj 



+ 



ki ■ k2 + (fci + fc2) -Pj 



m'jf -Bo(io) - - ml) (Bo(n) - 1) - -80(12) 



4a;^j Wife (Ci(io) - ^2(1) - 2£)oo(l)) 



Wfej Wife Bo (14) 



ki ■ Pj 



+ 2(fci • fca + {ki +/S2) -pj) 



+ 



fci • Pj 



mi 



mjwLLUik Bo{l) + 



2 (rnjj — 

ruj — m.f 



2 m^jj mj Lul^ ujik Bo{io) 



9 9 

m - — rnf 



+ 



rn| — m? 



( Hfe - rn]j) 



J "^kj 



u>ik - TOi Wfcj) Bo(ii)) - 2mj 



rnf + ml) lo^j ojik -80(12) + (m^ (2 mj mfc a;*fe 



+ (m| + - mlj) Wik) + {mi + ml) u-^ - m% m^ 0)*^. + 2 mfe Uik) Ukj) -80(13) j + 2 rrij 

X (m, uj*k [mk u^lj + mj ujkj) + i^ik {mk {mk wlj + mj ujkj) - mjj iVkj) ) Co(l) + 4 mj lv^j ujik Coo(l) 

+ 4mj a;feja;jfeCoo(3) - 2m| {mjwljUiik - miW*kU}kj) Ci{l) + 2mjmi {mjUi-kUJkj- {mi u^j uiik)) C2(l) 

+ \ TT, V ( i^i ^kj {{m-H -ml- ml) Uik - 2 m, mfc uj*k) - mj { {mf + ml) uj*k 

mi [mf -m^ -2ki-pj) y v j w 

-m%(j*k + 2mimkOJik)uJkj^ Bo{l) - {m'j -m"^) {m% - ml) u;*ki^kj Bo{ii) + mi (2mj {mf-mj) 

X u}*kj Wife -80(^2) + {rrij uj*kj (2 m, mfc uj*k + (m| + ml- mjj) uJik) + {mi {m] + ml) uJik - mjj mi uj*k 

+ 2m'^jmkL0ik) LOkj^ Baii^) - 2m^ (m,| - mf) (uj^k {m% luI^ Co(l) - mlcul^ Co(l) - m^ mtujkj Co(l) 
- 2ojIj Coo(l) + m'^jUjlj C'i(l)) + mj u^j C2(l) - a;*^. {mkcol^ Co(l) + mj Ukj (Co(l) + Ci(l) 



+ C2(1))))) 



+ 



mg {ml - ml) uj*^ujkj {Bo{l) - Bojii)) 
mi (m^ + 2ki- pj) 



{mjj - ml) ujl^LOjkBoih) 
m2 + 2 (fci • k2 + {ki + k2) ■ Pj) 



2oj*kjUjkm^ (Di2(l) + 1)23(1)), 

^Ij (■■-!) + 0(5) + r',.!---!) - 2 i}i2(i) /vi -pj) 



8 fci • Pj 

_l_ fcl • Pj fc2 • Pj 



fcl • Pj 

4 (mfc Wfcj Co(3) + mj w^j (Co(3) + Ci(5) + C2(5)) ) + 8 (mfc Wfcj Do(l) - mj ujI^ D2{1)) ki ■ pj 
mj {m]j - ml) uj^^ Bo{ii) + Bo{h) {mj ((m| + - m^) Wfcj + 2 mj mfc Wfcj) 



5(^1,^2) 

+ 2 (mj Wfcj + 2 mfc ojkj) {ki ■ k2 + {ki + k2) ■ Pj) ) 



im^ 
2 fci • ^2 



3 + 



2ki-k2 
2m'^ 



+ 6m^ 1 



2ki-k2 
2m'^ 



Co(i; 



^2 Ar 8mf 



(2 + (4m2-2/ei-/e2)Co(J/)). 



(A.1) 
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